Breit-Wigner resonances and the quasinormal modes of anti-de Sitter black holes 
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The purpose of this short communication is to show that the theory of Breit-Wigner resonances can be used 
as an efficient numerical tool to compute black hole quasinormal modes. For illustration we focus on the 
Schwarzschild anti-de Sitter (SAdS) spacetime. The resonance method is better suited to small SAdS black holes 
than the traditional series expansion method, allowing us to confirm that the damping timescale of small SAdS 
black holes for scalar and gravitational fields is proportional to r^^'~^, where is the horizon radius. The 
proportionality coefficients are in good agreement with analytic calculations. We also examine the eikonal limit 
of SAdS quasinormal modes, confirming quantitatively Festuccia and Liu's (3| prediction of the existence of 
very long-lived modes in asymptotically AdS spacetimes. Our results are particularly relevant for the AdS/CFT 
correspondence, since long-lived modes presumably dominate the decay timescale of the perturbations. 



I. INTRODUCTION 

It is well known that quasi-bound states manifest them- 
selves as poles in the scattering matrix, and as Breit-Wigner 
resonances in the scattering amplitude. Chandrasekhar and 
Ferrari made use of the form of the scattering cross sec- 
tion near these resonances in their study of gravitational-wave 
scattering by ultra-compact stars In geometrical units 

(c = G = 1), the Regge- Wheeler potential V{r) describing 
odd-parity perturbations of a Schwarzschild black hole (BH) 
of mass AI has a peak at r ~ 3M. Constant-density stellar 
models may have a radius R/M < 3 (but still larger than the 
Buchdal Hmit, R/M > 2.25). When R/M < 2.6, the ra- 
dial potential describing odd-parity perturbations of the star 
(which reduces to the Regge-Wheeler potential for r > R) 
displays a local minimum as well as a maximum. If this min- 
imum is sufficiently deep, quasi-stationary, "trapped" states 
can exist: gravitational waves can only leak out to infinity by 
"tunneling" through the potential barrier. Since the damping 
time of these modes is very long, Chandrasekhar and Ferrari 
dubbed them "slowly damped" modes iQl]. 

For trapped modes of ultra-compact stars the asymptotic 
wave amplitude at spatial infinity 'I' ~ a cos uir + f3 sin uur 
has a Breit-Wigner-type behavior close to the resonance: 

+ w const [(w - Wii)^ + w|] , (1.1) 

where ujJ^ is the lifetime of the quasi-bound state and ivj^ its 
characteristic "energy". The example of ultra-compact stars 
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shows that the search for weakly damped quasinormal modes 
(QNMs) corresponding to quasi-bound states {oj ~ uij^ — iojj 
with ujj ^ ojfj) is extremely simplified. We locate the reso- 
nances by looking for minima of + fP on the uj — oj^. line, 
and the corresponding damping time loj can then be obtained 
by a fit to a parabola around the minimum JllQl- 

Here we show that this "resonance method" can be used 
very successfully in BH spacetimes. The resonance method 
is particularly valuable in studies of asymptotically AdS BHs. 
The QNM spectrum of AdS BHs is related to thermalization 
timescales in a dual conformal field theory (CFT), according 
to the AdS/CFT conjecture ij^. Analytic studies of wave scat- 
tering in AdS BHs have previously hinted at the existence of 
resonances (see Fig. 9 in Ref. |01); here we show that these 
are indeed Breit-Wigner resonances. 

Various analytic calculations recently predicted the exis- 
tence of long-lived modes in asymptotically AdS BH space- 
times m 1^ 01- These modes will presumably dominate 
the BH's response to perturbations, hence the thermalization 
timescale in the dual CFT. Since their existence may be very 
relevant for the AdS/CFT conjecture, we decided to investi- 
gate numerically these long-lived modes. In Section we 
confirm the existence of quasi-bound states for small SAdS 
BHs, first predicted by Grain and Barrau U], partially correct- 
ing some of their predictions. In Section|III]we re-analyze the 
eikonal limit of SAdS QNMs studied by Festuccia and Liu 
10], finding excellent agreement with their calculations. 

It may be useful to point out that a different but intimately 
related method (the complex angular momentum approach, a 
close kin of the theory of Regge poles in quantum mechanics) 
has been used in the past to study QNMs in asymptotically 
flat BH spacetimes |8, 9:, Some aspects of the relation 
between the resonance method and the theory of Regge poles 
are illustrated, for example, in Ref. |01 ■ 
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II. QUASI-BOUND STATES IN SADS BLACK HOLES 

In this paper we will focus on SAdS BHs in four spacetime 
dimensions, but our results are trivially extended to higher di- 
mensions. Scalar (s — 0), electromagnetic (s = 1) and vector- 
type (or Regge-Wheeler) gravitational perturbations (s = 2) 
of SAdS BHs are governed by a second-order differential 
equation for a master variable v]/ ^ [Tl|] : 



— + // — 



Vl,s = f 



(1 



'2M 



2L2 



(2.1) 
(?.2) 



where / ^r^/L^ + l — ro/r. Lis the AdS radius, ro is related 
to the horizon radius through vq/L ~ (r^/L)^ + r+/L, 
and we assume that the perturbations depend on time as e^'"*. 
As usual, we define a "tortoise" coordinate by the relation 
dr/dr^ = / (so that — > — oo as r ^ r+) . 
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FIG. I: Potential for scalar field (s = 0) perturbations of a SAdS 
background with I — 0. Different lines refer to different values of 
r+ /L. A potential well develops for small BHs (r+ /L < 1). 



correspondence ifl^fTslfTiifTsll . With this caveat, most calcu- 
lations in the literature assume Dirichlet boundary conditions, 
so we choose to work with those. 

In general, a solution with the correct boundary conditions 
at infinity behaves near the horizon (r, —00) as 

* - Ai„e~''^''- +Aoute"^'^-' ~ acosLur^+Painujr^ , (2.5) 

with a = Aout + , (3 = i(^out — ^in)- For increased 
numerical accuracy, in our calculations we use a higher-order 
expansion of the form 

An{l+a{r-r+))e-'^^'+Aout{l+a*{r-r+))e'^-- , (2.6) 

with 



2l{l + 1) + {s^ ~ 1) ((.s^ - 6)rl/L^ - 2) 
(2r+/L)(l + Srl/L'^ - 2iLur+) 



(2.7) 



The problem is analogous to axial gravitational-wave scat- 
tering by compact stars, as long as we replace the "outgoing- 
wave boundary condition at infinity" in the stellar case by 
an "ingoing-wave boundary condition at the horizon" in the 
SAdS case (compare our Fig. [T]with Fig. 1 inRef. [Ql]). Quasi- 
bound states for the potential i2.2i should show up as Breit- 
Wigner resonances of the form dl.lt for real us. 




The potential for scalar-field perturbations of SAdS BHs is 
shown in Fig. [T] for different values of r+/L, ranging from 
"large" BHs with r+ /L ^ 10** to "small" BHs with r+ /L ^ 
10^^. Notice how a potential well of increasing depth and 
width develops in the small BH limit (r+ /L <C 1). 

Close to the horizon, where the potential V; ^ 0, we re- 
quire ingoing-wave boundary conditions: 

^ ~ e"*"""* , ^ -00 (r ^ r+) . (2.3) 

Near spatial infinity (r — » 00) the asymptotic behavior is 

«',s=o - Ar-^ + Br , 2 A/r + B . (2.4) 

Regular scalar-field perturbations should have 5 = 0, cor- 
responding to Dirichlet boundary conditions at infinity. The 
case for electromagnetic and gravitational perturbations is less 
clear, and there are indications that Robin boundary condi- 
tions may be more appropriate in the context of the AdS/CFT 



FIG. 2: A plot of + 0^ for scalar field SAdS perturbations with 
I — 0, r+ / L — 10~^. Resonances are seen when lor ~ 3 + 2n, 
i.e. close to the resonant frequencies of the pure AdS spacetime. In 
the inset we show the behavior near the minimum, which allows us 
to extract the decay time by a parabolic fit. 

Fig. |2] shows a typical plot of + /S^ as a function of 
ijJbL. The pronounced dips correspond to the location of a 
resonance, oj ^ ujr, and the inset shows a zoomed-in view of 
one such particular resonance. Once a minimum in + 0^ 
is located, the imaginary part w/ can be found by a parabolic 
fit of + 0^ to the Breit-Wigner expression ( 11.11 ). Alterna- 
tive (but equivalent) expressions are loj = —j3/a' = a/ (3' , 
where a prime stands for the derivative with respect to lob., 
evaluated at the minimum filial- We used these three differ- 
ent expressions to estimate numerical errors in the computed 
quasinormal frequencies. In the next subsections we briefly 
report on our results for 5 = and s — 2, respectively. 
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A. Scalar field perturbations 

The series solution method presented by Horowitz and 
Hubeny ^ was used by Konoplya in Ref. lfT6ll to compute 
quasinormal frequencies of small SAdS BHs. The series has 
very poor convergence properties for r+/L < 1, and QNM 
calculations in this regime take considerable computational 
time. As seen in Fig. [T] the potential for small SAdS BHs is 
able to sustain quasi-bound states, so we expect the resonance 
method to be well adapted to the study of small BHs. 



TABLE 1: The fundamental I = QNM frequencies for small SAdS 
BHs for selected values of r+ / L. The series solution data is taken 
from Ref. IT^ . In the table, Slor L = 3 — ujr L. 

series resonance 

L/r+ 5ujRL'^/r+ a^/L^/r^ SujRL^/r+ u}iL^/r\ 



12 


3.064 


9.533 


2.992 


9.662 


20 


2.922 


7.720 


2.912 


7.840 


30 


2.805 


6.660 


2.802 


6.714 


50 






2.700 


5.952 


100 






2.610 


5.471 


200 






2.580 


5.266 


500 






2.560 


5.158 


1000 






2.554 


5.125 


2000 






2.550 


5.109 


5000 






2.550 


5.100 



In TableUwe list QNMs for / = scalar field perturbations 
and for different BH sizes, comparing (where possible) results 
from the resonance method with Konoplya's series expansion 
calculation. A cubic fit of our data for L/r_(. > 30 yields 
Luj L = 5.00r^/L2 + 47.70r^/L^, a quartic fit yields ajj = 
5.09(r+/L)2 + 33.59(r+/L)3 + 485.09(r+/L)4, and fits with 
higher order terms basically leave a and b unchanged with 
respect to the quartic fit. The numerical results are consistent 



with Horowitz and Hubeny's prediction that uji 



and 



they are in very good agreement with analytic predictions by 
Cardoso and Dias UtIi . who derived a general expression for 
the resonant frequencies of small BHs in AdS for Muiji <^ 1 
regime. Setting a = in Eq. (33) of Ref. UtIi . their result is 



Luj = I + 3 + 2n — iujjL . 
where n is a non-negative integer and 



(2.8) 



ui L 



-70 



(/ + 3 + 2n) (r+/L)''+' /ir 



(2.9) 



70 



2-1-61 (^1)2 p _ ^/2f r [5 + 2/ + 2n] 
(3 + 2n)(3 + 21 + 2n)r [I + l/2f T [2 + 2n] 



For Z = n = one gets uji = 16{r+/L)^/TT ~ 5.09{r+/Lf, 
in excellent agreement with the fits. For general / Eq. ( 12.8b 



predicts an 



„2+2l 



dependence, in agreement with our results 



for / = 0, 1, 2. Moreover we find excellent agreement with 
the proportionality coefficient predicted by Eq. i2.9i . Higher 
overtones are also well described by Eqs. (I2.8l l and ( 12. 9t . 

Our results show that the resonant frequency ujj^ always ap- 
proaches the pure AdS value in the small BH limit, generally 
confirming the analysis by Grain and Barrau |5j]. However our 



numerics disagree with Grain and Barrau's semi-classical cal- 
culation of the monopole mode (/ = 0). We find that all modes 
including the monopole reduce to pure AdS in the small BH 
limit. More precisely, as r+ /i ^ we find 

ujRL = l + 3 + 2n - kin r+/L , n = 0, 1, 2, . . . (2.10) 

with kio ~ 2.6, 1.7, 1.3 for / = 0, 1, 2 , respectively. 




FIG. 3: Track traced in the complex plane (ujrL,ujiL) by the fun- 
damental I = 0,1 scalar field QNM frequencies as we vary the BH 
size r+ /L. Counterclockwise along these tracks we mark by circles 
and diamonds the frequencies corresponding to decreasing decades 
in r+/L (r+/L = 10^ 10\ lO", 10"\ ...). 



Our results for scalar field perturbations are visually sum- 
marized in Fig. [51 where we combine results from the res- 
onance method and from the series solution (compare e.g. 
Fig. 1 of Ref. ifisll ). Modes with different I coalesce in the 
large BH regime (top-right in the plot), as long as I <C r+ /L. 
As shown in Fig.[Tl the potential for small SAdS BHs devel- 
ops a well capable of sustaining quasi-stationary, long-lived 
modes. It should not be surprising that small and large BH 
QNMs have such a qualitatively different behavior 



B. Gravitational perturbations 

We have also searched for the modes of Regge- Wheeler or 
vector-type gravitational perturbations. The potential ( 12.21 ) for 
gravitational perturbations of small BHs does not develop a 
local minimum. Nevertheless it develops a local maximum 
which, when imposing Dirichlet boundary conditions [see the 
discussion following Eq. ( 12. 41 )1 can sustain quasi-bound states. 
For this reason we expect the resonance method to be useful 
also in this case. 

For small BHs our numerical results agree with the fol- 
lowing analytic estimate, derived under the assumption that 
Mu;r < 1 tl9]: 



w/L~ -72(/ + 2 + 2?i) (r+/L) 



21+2 



(2.11) 



with 
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{I + 1)(2 + / + 2n) r [1/2 - /] r [if T[3 + l]r [2(2 + l + n)] 
{I - l)2(2n)!r [Z + 3/2] r [2Z + 1] r [2^ + 2] r [-1/2 - Z - 2n] r [7/2 + ; + 2n] 



(2.12) 



r 



For 71 = and / = 2 this implies ivj = 1024(r+/L)6/457r ~ 
7.24(r+/L)®, while a fit of the numerics yields ojj ~ 
7.44(r+/L)^. Again, ujuL approaches the pure AdS value 
in the limit / L ^ 0: 

iURL^l + 2 + 2n- kin r+lL , n = 0, 1, 2, . . . (2.13) 

For the fundamental I = 2 mode we find fc2o ^ 1.4. 

III. LONG-LIVED MODES IN THE EIKONAL LIMIT 

A recent study of the eikonal limit (/ > 1) of SAdS BHs 
suggests that very long-lived modes should exist in this regime 
10]. Define rt > Tc to be the two real zeros (turning points) of 
cj|j. — f /r^ = 0. Then the real part of a class of long-lived 
modes in four spacetime dimensions is given by the WKB 
condition 



The resonance method is well suited to analyze the eikonal 
limit, especially for small BHs (for large BHs the existence 
of a dip in the potential well requires very large values of 
I, which are numerically hard to deal with). In Table HIl we 
compare the WKB results against numerical results from the 
resonance method. Unfortunately machine precision limita- 
tions do not allow us to extract extremely small imaginar 
parts. The agreement with the WKB condition of Ref 
is remarkable, even for relatively small values of /. Our nu- 
merics conclusively confirm the existence of very long-lived 
modes in the SAdS geometry, but the numerical results for the 
damping timescales disagree by orders of magnitude with the 
corresponding results by Grain and Barrau A reanalysis 
of the assumptions implicit in their method would be useful to 
understand the cause of this disagreement. 



ary 

i 



P f 



dr = 7r(2n + 5/2) , 



rf 

' I'b 'J 

where p = I + 1/2. Their imaginary part is given by 



(3.1) 



8ll>r 



logT = 2i 



rf 



dr . (3.2) 



The pref actor 7, not shown in Ref. i^, can be obtained by 
standard methods 1I20I uJJl with the result 



cos-2 X 



dr 



■ dr 



(3.3) 



IV. CONCLUSIONS AND OUTLOOK 

The method described here provides a reliable and accu- 
rate alternative to the series solution method |3J, to be used 
in regimes where the former has poor convergence properties. 
Together, the two methods allow an almost complete charac- 
terization of the spectrum of BHs in AdS backgrounds, en- 
compassing both small and large BHs. As an application of 
the method, we have explicitly confirmed the existence of the 
weakly damped modes predicted by Refs. IHI^. 

Extensions of the resonance method to higher-dimensional 
^ and charged geometries ifisi I22I1 should be trivial. Our 
techniques may be useful to verify the existence of the highly- 
real modes predicted by Ref. |01. Finally, it would be inter- 
esting to investigate whether the resonance method described 
here is of any use to investigate the eikonal limit of QNMs in 
asymptotically flat BH spacetimes. 



TABLE II: The QNM frequencies for a r+/L = 0.1 BH for selected 
values of I, and comparison with the corresponding WKB prediction. 



{l,n) 



WKB 



resonance 

ujrL logoj/L 



3,0 


5.8668 


-16.40 


5.8734 


-17.03 


3,1 


7.6727 


-12.61 


7.6776 


-12.97 


3,2 


9.4189 


-9.40 


9.4219 


-9.60 


4,0 


6.8830 


-22.13 


6.8889 


-22.84 


4,1 


8.7139 


-18.30 


8.7184 


-18.41 


4,2 


10.4960 


-14.65 


10.4996 


-14.72 


5,0 


7.8945 


-28.64 


7.8997 


-28.76 


5,1 


9.7426 


-23.97 


9.7466 


-24.02 


5,2 


11.5482 


-20.03 


11.5516 


-20.06 
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